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LONG-TERM  GOALS 

The  goal  of  this  research  is  to  understand  the  effects  of  range-dependent  sediment  properties  on  the 
acoustic  field  in  2-D  shallow  water  environments.  This  information,  in  part,  is  required  to  solve  the 
statistical  inference  problem  in  inhomogeneous  shallow  water  environments. 

OBJECTIVES 

The  objective  of  the  current  work  is  to  create  and  benchmark  a  two-way  coupled  mode  model  that  will 
later  be  used  to  analyze  the  effects  of  range-dependent  sediment  properties  on  the  acoustic  field  in 
typical  littoral  environments.  An  accurate  coupled  mode  model  is  required  to  properly  capture  the 
physics  of  shallow  water  acoustic  propagation  and  to  classify  the  propagation  as  belonging  to  either  the 
adiabatic  or  coupled  mode  regime.  It  is  of  particular  interest  to  classify  acoustic  data  taken  during  the 
Shallow  Water  2006  (SW06)  experiment  as  being  adiabatic  or  coupled  mode  and  to  identify  the 
physical  mechanisms  that  cause  mode  coupling. 

APPROACH 

The  model  developed  in  this  work  is  based  on  two-way  coupled  mode  theory.  The  inhomogeneous 
Helmholtz  equation  for  pressure  pa(z,r)  due  to  a  continuous  wave  (CW)  point  source  located  at  (zs-  rs) 
is 


p(z,r)V- 


1 

_p(z,r ) 


vPco(z,r) 


+  k2(z,r)pco(z,r ) 


-471 Po  8(r-rs)d(z-Zs), 


(1) 


where  k  =  2nf/c{z,r)  is  the  spatially-dependent  acoustic  wavenumber,  z  is  the  depth  coordinate,  r  is 
the  horizontal  coordinate,  p  is  the  spatially-dependent  fluid  density,  and  po  is  the  CW  source  spectral 
amplitude  referenced  to  1  m.  The  fluid  properties  for  c  are  allowed  to  be  completely  arbitrary  in  depth 
and  range  while  the  properties  for  p  are  allowed  to  be  piecewise  continuous  in  depth  and  arbitrary  in 
range.  Continuity  of  pressure  at  interfaces  and  continuity  of  particle  velocity  normal  to  fluid  interfaces 
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are  properly  enforced  to  conserve  energy.1,2  The  solution  for  pco(z,r)  can  be  represented  by  the  normal 
mode  summation 

M 

Pco(z,r)  =  £  Rm(r\zs,rs)^m(z,f),  (2) 

m=  1 

where  M  is  the  total  number  of  modes  included  in  the  summation,  Rm(r\zs,  rs)  are  range-dependent 
modal  amplitudes  specific  to  the  source  position  zs  and  rs,  and  <j>m(z,r)  are  the  depth  dependent  local 
mode  functions.  Note  that  in  this  work,  the  possibility  of  horizontal  (out  of  plane)  refraction  is 
neglected.  Therefore,  the  range  coordinate  r  becomes  v  or  r  for  Cartesian  or  cylindrical  geometries, 
respectively. 

The  major  challenge  of  the  two-way  coupled  mode  model  used  in  this  research  is  to  solve  for  the 
range-dependent  modal  amplitudes  Rm(r\zs,  rs)  that  appear  in  Eq.  (2).  This  is  done  by  transforming  Eq. 
(1)  into  an  integral  equation  of  the  form 


Rm{r\zs,rs)  =  -4npo G,„ (r,  rs)  -  J Gm{r,r) 


M 


T  (Bmw(r/)  •  Vp  +  Cff 


r)Rn(r 


1^,5) 


n=  1 


dr' ,  (3) 


where  Gm(r,rs )  is  a  two-way  adiabatic  Green’s  function  and  Bmn(r)  and  Cmn  are  mode  coupling 
matrices.  The  solution  method  for  Eq.  (3)  is  relatively  efficient  when  using  the  Lanczos  routine 
proposed  by  Knobles.3  It  should  be  noted  that  the  coupled  mode  model  developed  here  is  capable  of 
including  both  the  trapped  and  continuum  components  of  the  modal  spectrum. 


WORK  COMPLETED 

Much  of  the  work  completed  was  to  implement  the  coupled  mode  model  in  MATLAB  and  then  to 
benchmark  the  new  model  against  other  known  models  and  accepted  results.  The  new  model  was 
benchmarked  for  a  number  of  cases,  including  problems  with  range-dependencies  in  the  bathymetry,  the 
water  column,  and  the  sediment  geoacoustic  properties.  Through  the  benchmarking  process,  the  model 
was  shown  to  be  in  agreement  with  other  model  solutions,  specifically  a  PE  model  and  a  two-way 
coupled  mode  model  based  on  the  differential  equation  form  of  the  2-D  Helmholtz  equation. 

The  integral  equation  coupled  mode  model  was  then  applied  to  a  typical  propagation  environment 
encountered  during  the  SW06  experiment  to  see  if  coupled  mode  effects  due  to  seabed  inhomogeneities 
were  present.  The  particular  environment  under  analysis  included  variable  bathymetry  and  sub-bottom 
interface  depths.  Range-independent  sediment  properties  and  water  column  sound  speed  profiles  were 
assumed  for  the  initial  analysis,  which  is  an  assumption  that  may  be  suspect  for  the  SW06  environment 
(consequently  the  initial  analysis  may  underestimate  the  amount  of  mode  coupling).  However,  even  in 
this  simplified  environment,  it  was  observed  through  the  model  calculation  that  coupled  mode  effects 
may  be  present  in  the  measured  acoustic  data. 

An  initial  attempt  was  made  to  match  the  modeled  data  with  measured  acoustic  data  on  one  of  the 
vertical  line  arrays  (VLA)  deployed  during  SW06.  A  good  initial  fit  between  model  and  data  was 
sought  in  order  to  make  conclusive  statements  about  whether  coupled  mode  effects  are  present  in  the 
measured  acoustic  data.  The  effect  of  horizontal  variability  in  the  environment  was  shown  to  have  a 
minimal  effect  on  the  differences  between  modeled  and  measured  transmission  loss. 
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RESULTS 


This  section  presents  results  for  (A)  the  comparison  of  the  coupled  mode  model  to  other  propagation 
models  for  selected  benchmark  problems  and  (B)  model-data  comparisons  for  a  propagation  track  in  the 
SW06  experiment. 

A.  Benchmark  comparison 

The  results  of  two  benchmark  problems  are  shown  here.  Throughout  this  section  the  new  integral 
equation  coupled  mode  model  will  be  designated  by  the  acronym  IECM,  Collin’s  range-dependent 
acoustic  model4  will  be  designated  by  RAM,  and  Evan’s  coupled  mode  model5  will  be  designated  by 
COUPLE.  Comparisons  to  RAM  and  COUPLE  were  made  as  these  programs  represent  the  current  state 
of  the  art  in  range  dependent  propagation  modeling. 

1.  The  shelf  break  front  problem  was  proposed  as  part  of  the  SWAM99  workshop.  The  benchmarked 
case  shown  here  is  after  the  slightly  modified  problem  published  by  Knobles.6  Both  the  bathymetry  and 
water  column  sound  speed  profile  are  range-dependent.  The  comparison  between  the  IECM  and  RAM 
solutions  are  shown  in  Fig.  1,  with  contour  lines  drawn  at  3  dB  intervals  throughout  the  range-depth 
plane.  There  is  excellent  agreement  between  the  two  models  in  the  water  column  but  discrepancies  exist 
in  the  field  inside  the  sediment  layers.  As  is  typical  when  using  the  Pekeris  branch  cut  (which  is  done  in 
this  research  when  using  ORCA  to  find  the  local  mode  depth  functions),  the  unphysical  behavior  of  the 
leaky  modes  in  the  sediment  layers  cause  the  field  in  the  sediment  to  be  in  error. 

2.  Outing7  proposed  a  benchmark  problem  involving  a  downslope.  His  Example  A  was  slightly 
modified  for  the  benchmark  case  shown  here.  A  source  is  located  at  zs  =  10  m  in  the  water  column 
emitting  a  CW  sound  at  /  =  100  Hz.  The  comparison  between  the  IECM,  RAM,  and  COUPLE 
solutions  are  shown  to  be  in  excellent  agreement  in  Fig.  2. 


Figure  1:  Top:  Contour  plot  of  TL  for  the  shelf  Figure  2:  Top:  Contour  plot  of  TL  for  Outing’s 
break  front  problem.  Bottom:  TL  at  zr  =  35  m.  Example  A  problem.  Bottom:  TL  at  zr  =  20  m. 


The  IECM  model  was  applied  to  an  environment  that  occurred  during  the  SW06  experiment  as  the  R/V 
Knorr  towed  a  J-15  source  from  39.0457  N,  73.1362  W  to  38.9923  N,  73.0172  W  on  August  26th 
between  the  hours  of  02:15  to  03:30.  The  bathymetry  underneath  the  Knorr  was  measured  by  an 
onboard  Knudsen  echosounder  and  is  shown  in  Fig.  3.  The  variation  in  the  bathymetry  over  the  course 
of  the  track  was  approximately  20%  of  the  mean  water  depth.  The  measured  bathymetry  was  used  to 
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create  a  representative  acoustic  propagation  track  between  source  and  receiver.  Nominal  sound  speed, 
density,  and  attenuation  values  of  1650  m/s,  1.6  g/cm3,  and  0.08  dB/A  and  1725  m/s,  2.0  g/cm3,  and 
0.05  dB/A  were  used  for  the  first  layer  and  halfspace,  respectively.  The  average  water  column  sound 
speed  profile  measured  on  the  nearby  Shark  array  over  the  time  of  the  tow  was  used  in  the  model. 


* -  Bathymetry  measured  by  the  onboard  Knudsen  12  kHz  echosounder 
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Figure  3:  Bathymetry  as  measured  by  the  onboard  Knudsen  echosounder  for  a  down-shelf  track 

from  SWAMI 32. 

Transmission  loss  as  a  function  of  range  was  computed  for  three  different  frequencies  (53,  103,  and  203 
Hz)  and  three  different  receiver  depths  (10.49,  28.36,  and  46.22  m).  Both  the  adiabatic  and  coupled 
mode  solutions  were  computed  with  the  IECM  model.  The  comparisons  of  these  two  solutions  are 
shown  in  Fig.  4.  It  is  observed  that  the  bathymetry  variation  present  in  this  environment  is  sufficient  to 
cause  coupled  mode  propagation  at  these  frequencies.  It  is  also  observed  that  the  degree  of  mode 
coupling  increases  with  frequency. 


Figure  4:  Adiabatic  (red)  and  IECM  (blue)  TL  solutions  for  the  SW06  environment  at  53, 103,  and 

203  Hz  and  three  different  receiver  depths. 


B.  Model-data  comparison 

Through  the  efforts  described  above,  it  was  discovered  that  coupled  mode  effects  may  be  present,  and 
even  significant,  in  the  SW06  environment.  An  effort  was  made  to  observe  these  effects  in  measured 
acoustic  data.  This  first  required  a  more  careful  construction  of  the  acoustic  propagation  path  between 
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source  and  receiver  in  order  to  decide  whether  misfit  in  the  model-data  transmission  loss  was  caused  by 
horizontal  inhomogeneities.  As  is  commonly  done  in  propagation  modeling,  a  single  acoustic  track  was 
constructed  between  the  receiver  location  (the  SWAMI32  VLA  located  at  39.0611  N  73.1308  W)  and 
the  final  position  of  the  source  along  the  tow  track  (38.9923  N,  73.0172  W).  The  common  assumption  is 
that  horizontal  variability  in  the  environment  is  small  enough  to  be  neglected.  This  track  will  be  called 
the  “endpoint”  track  and  is  shown  in  blue  in  Fig.  5.  The  actual  acoustic  tracks  propagate  from  each 
source  position  on  the  red  track  back  to  the  SWAMI32  VLA.  Each  separate  acoustic  track  is  plotted  in 
red  on  the  bottom  of  Fig.  5  to  illustrate  horizontal  bathymetry  variations  in  this  region. 


Figure  5:  Top:  Knorr  ship  track  (red)  and  endpoint  acoustic  track  (blue)  plotted  on  top  of  bathymetry 
contours  (grayscale).  CTD  measurements  shown  as  dots.  Bottom:  Bathymetry  and  sub-bottom  layer 
depths  for  multiple  acoustic  tracks  (red)  and  endpoint  acoustic  track  (blue). 

Because  of  its  computational  efficiency,  the  propagation  model  RAM  was  used  to  compute  the 
transmission  loss  along  1 80  different  acoustic  propagation  paths  between  the  receiver  and  the  moving 
source.  RAM  was  also  used  to  compute  the  transmission  loss  along  the  single  endpoint  track.  The 
comparison  between  measured  and  modeled  transmission  loss  at  103  Hz,  while  using  nominal  values 
for  the  sediment  geoacoustic  parameters,  is  shown  in  Fig.  6.  The  degree  of  misfit  between  model  and 
data  does  not  change  significantly  at  103  Hz  by  considering  the  horizontal  variability  in  the 
environment.  This  is  an  important  result  for  this  particular  track  and  suggests  that  the  horizontal 
variability  is  small  enough  to  neglect  at  and  below  this  frequency.  This  will  allow  for  a  single-track 
computation  with  the  IECM  model  (which  is  much  more  computationally  intensive  than  RAM)  to 
analyze  the  significance  of  the  mode  coupling.  However,  it  is  possible  that  the  effect  of  horizontal 
variability  will  increase  with  frequency.  This  is  an  important  consideration  for  future  work,  given  that 
the  degree  of  mode  coupling  was  also  observed  to  increase  with  frequency. 
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Figure  6:  Comparison  between  measured  data  (black)  and  modeled  data  using  180  propagation 
paths  (red)  and  the  endpoint  path  (blue)  at  five  receiver  depths  on  the  SWAMI32  VLA. 

IMPACT/APPLICATIONS 

The  impact  of  this  research  is  to  increase  the  understanding  of  physical  propagation  mechanisms  in 
shallow  water  environments  and  to  identify  important  modeling  considerations  that  must  accounted  for 
to  decrease  uncertainty  in  model  prediction.  The  research  can  be  applied  to  propagation  models  of  all 
types  and  can  also  be  used  in  geoacoustic  inversion  and  inference  problems. 

TRANSITIONS 

The  primary  transition  for  this  project  is  an  accurate  2-D  coupled  mode  model  for  acoustic  propagation 
in  shallow  water  environments,  and  a  methodology  for  statistical  inference  for  seabed  parameters  in 
environments  that  have  horizontal  and  range  variability. 

RELATED  PROJECTS 

None. 
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